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Abstract 

In this paper we prove exact forms of large deviations for local times and intersection local 
times of fractional Brownian motions and Riemann-Liouville processes. We also show that 
a fractional Brownian motion and the related Riemann-Liouville process behave like constant 
multiples of each other with regard to large deviations for their local and intersection local times. 
As a consequence of our large deviation estimates, we derive laws of iterated logarithm for the 
corresponding local times. The key points of our methods: (1) logarithmic superadditivity of a 
normalized sequence of moments of exponentially randomized local time of a fractional Brownian 
motion; (2) logarithmic subadditivity of a normalized sequence of moments of exponentially 
randomized intersection local time of Riemann-Liouville processes; (3) comparison of local and 
intersection local times based on embedding of a part of a fractional Brownian motion into the 
reproducing kernel Hilbert space of the Riemann-Liouville process. 
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1 Introduction 



Let B H (t), t > be a standard (i-dimensional fractional Brownian motion with index H G (0,1). 
That is, B H (t) is a zero-mean Gaussian process with stationary increments and covariance function 



E 



B H (t)B H (s) T ] = ±{\t\ 2H + \s\ 2H -\t- s\ 2H } I d 



where 1^ is the identity matrix of size d. B H (t) is also a self-similar process with index H. The 
local time Lf(B H ) of B H (t) at x £ M. d is defined heuristically as 

Lf(B H )= f 5 x (B H {s)) ds, t>0. 
Jo 

It is known that Lf(B H ) exists and is jointly continuous in (t, x) as long as Hd < 1. By the 
self-similarity of a fractional Brownian motion, Lf(B H ) = t 1 ^ Hd L^ t (B H ). In particular, 

L° t (B H ) I t l - Hd L\(B H ). (1.1) 

Our first goal is to investigate large deviations associated with tail probabilities of L^(B H ). By the 
scaling given above, we may consider only t = 1. In the classical case, when H = 1/2 and d = 1, 
it is well known, see the book of Revuz and Yor 03], p240, that L^(B 1 / 2 ) = \U\ with U ~ M(0, 1). 
Consequently, 

lim a^logPfL^S 1 / 2 ) > a\ = --. 

a— >oo I J 2 

In Theorem 12.11 we prove that for a fractional Brownian motion a nontrivial limit 

lim a~ 1 / Hd \ogF{L° 1 {B H ) > a} 

a— >oo 

exists and we give bounds for this limit. 

Closely related to the fractional Brownian motion is the Riemann-Liouville process W H (t) with 
index H > which is defined as a stochastic convolution 

W H (t)= f\t-s) H - 1/2 dB(s), t>0, (1.2) 
Jo 

where B(t) is a d-dimensional standard Brownian motion. {W H (t)}t>o is a self-similar zero-mean 
Gaussian process with index H, as is B H (t), but W H (t) does not have stationary increments and 
there is no upper bound restriction on index H > 0. If L° (W ) denotes the local time of W H (t) 
at 0, then by the self-similarity we also have 

L° t (W H ) = t l - Hd L\{W H ). (1.3) 

The relation between W H (t) and B H (t) becomes transparent when we write a moving average 
representation of B H (t), t G M, in the form 

B H (t) = c H ! \{t - s) 11 - 1 ' 2 - (- S )+" 1/2 j dB(s), (1.4) 
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where 

c H = \/2H2 H B(l-H,H + l/2)~ 1/2 , (1.5) 

and B(-, •) denotes the beta function. The analytic derivation of ch is given for completeness in the 
Appendix (a different but equivalent form of ch is also derived in Mishura, |38[ Lemma A.0.1], by 
a Fourier analytic method). From ([1.4p we have a decomposition 

c H 1 B H (t) = W H (t) + Z H (t), (1.6) 

where 

Z H (t) = f \{t - s) H - 1 / 2 - {sf- 1 ' 2 ] dB{s) (1.7) 
is a process independent of W (t). 

This moving average representation for fractional Brownian motion was introduced in the pioneering 
work of Mandelbrot and Van Ness [36j and used extensively by many authors, sometimes with 
different normalizing constant ch in (jl.5p (e.g., Li and Linde |32] uses T(H + 1/2) -1 for c#). 

We will show that paths of Z H (t), away from t = 0, can be matched with functions in the reproducing 
kernel Hilbert space of W (t) (Proposition 13.51 Section [3,2p , This and the independence of Z H (t) 
from W (t) will allow us to show that large deviation constants of tail probabilities of L\(W H ) and 
of LKcJj 1 B H ) = c^jLl{B H ) are the same (Theorem 12. 2p . In this context we also want to mention 
Theorem 3.22 of Xiao, [47], who established bounds for tail probabilities of the local time L\ of the 
general Gaussian processes in the form 

-d < liminf — !— log{L? > a} < limsup — [ — logjL? > a} < -C 2 

a^co (p[a) a-^oo 4>{a) 

and raised a question on the existence of the limit (Question 3.25, |47|). Further, we cite the paper by 
Baraka, Mountford and Xiao (|4J) for some similar tail estimate of the local time of multi-parameter 
fractional Brownian motions. 

Next we will consider p independent copies B^ (t), . . . , Bp (t) of a standard d-dimensional fractional 
Brownian motion B H (t). Throughout this paper 

P* -=P/{P- 1) 

will stand for the conjugate to p > 1. Our next and main goal is to investigate large deviations 
for intersection local time a H (•) of B^(t), • • • , Bp (t), which is a random measure on (M + ) p given 
heuristically by 

a H (A) = / I J 8 (B?{ Sj ) - Bf +l {s j+l )) d Sl --- ds p , A c (M + ) p . 

J A j=1 

Quantities measuring the amount of self-intersection of a random walk, or of mutual intersection 
of several independent random walks, have been studied intensively for more than twenty years, 
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see e.g. |17| . [30] . |29| . |37| . |23j . |10) . [11] . This research is motivated by the role these quantities 
play in quantum field theory, see e.g. [18j, in our understanding of self-avoiding walks and polymer 
models, see e.g. [35], [25] . or in the analysis of stochastic processes in random environments, see 
e.g. [24J [20], [2], [19] , In the latter models dependence between a moving particle and a random 
environment frequently comes from the particle's ability to revisit sites with an attractive (in some 
sense) environment. Consequently, measures of self-intersection quantify the degree of dependence 
between movement and environment. Typically, in high dimensions, this dependence gets weaker, 
as the movements become more transient and self-intersections less likely. Investigation of large 
deviations for intersection local times is closely related to asymptotics of the partition functions in 
above models. 

There are two equivalent ways to construct a H (A) rigorously. In the first way, a H (A) is defined as 
the local time at zero of the multi-parameter process 

X(h,--- ,t p ) = (flf (ii) - (t 2 ), • • • ,B^_ x {t p ^)-B^{t p )) (*!,-•■ ,t p ) e (R+r (1.8) 

More precisely, consider the occupation measure 

HA{B)= [ 1 B {B?(8 X ) - B?(s 2 ),- •• ,BjH 1 (s p - 1 )-Bj?(s p ))ds 1 ---ds p , Bel*' 1 *. 
J A 

By Theorem 17. 1( as Hd < p* , there is a density function a H (A, •) of ^a(') such that if A = 
[0,h] x ••• x [0,t p ], then a H ([0,h] X • • • X [0, t p ], x) is jointly continuous in (tx, ■ ■ ■ , t p , x). We define 
a H (A) :=a H (A, 0). 

For the second way of constructing a H (A), write for any e > 



af (A) = / / \[p e (Bf ( Sj )-x)dsi---ds p dx, (1.9) 

JR d J A jj[ 

where p e are probability densities approximating 5q as e — > 0. Notice that 



af(A)= f h e (B*( Sl )-B?(s 2 ),--- ,B*_ 1 (s p ^ 1 )-B p H (s p ))ds 1 ---ds 
J A 

h e (x)a H (A, x)dx 



where 

„ p-i p-1 

h e (xi,--- = / Pe(-X) T\pA ~ X 

3=1 k=j 

is an probability density on ^ approaching (5q(x^, • • • ^x p —\ ) as e — > + . 

By the continuity of a H (A, x), lim e _ s> o+ a e^(^) = ct H (A) almost surely. Applying Proposition 13.11 
to the Gaussian field given in (|1.8|) . the convergence is also in C m for all positive m. This way of 
constructing a H (A) justifies the symbolic notation 

a H (A)= f [ T\S Q (Bf ( Sj ) -x) dsf-dspdx. 

jR d J A j =1 
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In the special case p = 2 and Hd < 2. Nualart and Ortiz-Latorre |39| proved that a^QO, t{\ x [0, £2]) 
converges in C 2 as e — )• + , with 

p e (x) = (2evr)- d/2 exp{-|x| 2 /2e}. (1.10) 
For the Riemann-Liouville process W H (t) an analogous construction of the intersection local time 




can be done under the same condition Hd < p* . 

By the self-similarity of B H (t) and W (t), for any t > 

a H ([0,t} p ) =t p - Hd( P-^a H ([0,l] p ) (1.11) 

and 

a H ([0,t] p ) =t p - Hd ( p -Va H ([0,l} p ). (1.12) 

Finally, we would like to discuss this research in a more general context of Markovian versus non- 
Mar kovian structures. Naturally, most of the existing results on large deviation for (intersection) 
local time have been obtained for Markov processes such as Brownian motions, Levy stable processes, 
general Levy processes, and random walks. The underlying Markovian structure has been essential 
for the methods in these studies; see Chen |llj for references and a systematical account of such 
works. Departures from Markovian models are often driven by the underlying physics to match 
the required level of dependence (memory) and smoothness/roughness of sample paths. Fractional 
Brownian motion and Riemann-Liouville processes are the most natural candidates as extensions of 
Brownian motion into the non-Markovian world. They offer the existence of the intersection local 
time for any number p of processes in any dimension d as long as H is sufficiently small. Therefore, 
they may help scientists to build more realistic and robust models while posing serious challenge to 
mathematicians due to the non-Markovian nature. 

In this paper, we mainly use Gaussian techniques motivated from the study of continuity properties 
of local time, and more generally, from theory of Gaussian processes. It is also helpful to see con- 
nections between small ball probability estimates and tail behavior of the local time. Indeed, large 
value of the local time at zero means that the process stayed for a long time in a small neighborhood 
of zero. By this analogy, Propositions 13.11 and 13.31 can be motivated by the corresponding results 
for small balls (see comments preceding these propositions in Section \3.1\i . 
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2 Main results 



Theorem 2.1 Let B (t) be a standard d- dimensional fractional Brownian motion with index H 
such that Hd < 1. Then the limit 

lim aT 1 ^ log ¥{L\{B H ) >a} = -9{H, d) (2.1) 

a— >oo 

exists and 9(H, d) satisfies the following bounds 

{txc 2 h /H) 1I(2H) 6 (Hd) < 9(H, d) < (2n) 1 ^ 2H ^e (Hd), (2.2) 
where ch is given by U.5\) and 

*w=«(W^r) • (2 - 3) 

Notice that in the classical case of one-dimensional Brownian motion, (|2.2p becomes the equality. 
The fact that the lower bound is less than or equal to the upper bound in (I2.2j) is equivalent to 
cL < 2H, which can also be seen directly. Indeed, from (3.16) 

% = Vax(B H {l)\B H (s), s < 0) < V&r(B H (l)) = 1. (2.4) 
The equality only holds for a Brownian motion, i.e., H = 1/2. 

Theorem 2.2 Let W H (t) be a d-dimensional Riemann-Liouville process as in (11. 2p such that Hd < 
1. Then the limit 

lim a -1/(H * \og¥{L\(W H ) >a} = -9(H,d), (2.5) 

a— >oo 

exists with 

6(H,d) = (c H y 1/H 6(H,d), (2.6) 
where 9(H,d) is as in Theorem \2.1\ and ch is given by U.5\) . 

Theorem 2.3 Leta H (-) be the intersection local time of p-independent d-dimensional Riemann-Liouville 
process W^{t), • • • , W^(t), where Hd < p*. Then the limit 

lim a- p */^ dp )logP{a^([0,l] p ) > a) = -K(H,d,p) (2.7) 
exists and K(H,d,p) satisfies the following bounds 

p— (l-_ ) Hd [-) 2 V^r(l-— ) Hd <K(H,d,p) (2.8) 



where ch is given by U.5\) . 



V 

Hd 
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There is a direct way to show that the lower bound is less than or equal to the upper bound in 
flZU) . Observe that by Holder inequality 1 + t 2H > p l /P(jp*) l IP* t 2H l p * which leads to 

poo 

/ (i + t^Y^e-Ut < p-"/(2p)( p *)-«'/(2p*) r ( 1 _ Hd/p *y 
Jo 

After cancellation on both sides of (|2.8p . the problem is then reduced to examining the relation 
c% < 2H, which is given in (12. 4ft , 

Theorem 2.4 Let a H (-) be the intersection local time of p '-independent standard d- dimensional 
fractional Brownian motions B^(t), ■ ■ ■ ,B^{t), where Hd < p*. Then the limit 

lim a' pt ^ Hdp hog¥\a H ([0,l] p ) > «} = -K(H,d,p) (2.9) 

exists with 

K(H,d,p) = c]i H K(H,d,p). (2.10) 

Our results seem to be closely related to the large deviations of the self-intersection local times 
heuristically written as 

P H {M P <) = / HSo{B H (s j )-B H (s j+l ))ds 1 ---d Sp 

J \o,t] p < j= i 

where 

[0,i]* = {(si,-"" ,s P ) € [0,t] p ; Sl < ••• < s p }. 
In the case when Hd < 1, we can rewrite 



/^(M P <)= -J [LUB H )] P dx. 



1/p 



\ VP 



To see the connection between a H and /3 , notice that by Holder inequality and arithmetic and 
geometric mean inequality 

l II ^{Bf)dA <\t(l d [Lt(Bf)Yd X ) 

Thus, for any 6 > 

Eexp{^a £ W ! ^([0,l] p )) 1/P } < Eexp ^Qp~ 1 a E ~ B ^rT f J [Lf(B H )] p dx^j 

On the other hand, by Theorem 12.41 and Varadhan's integral lemma, 

Hm a - p */( Hd rt log Eexp {Op^a*^ (a H '([0, 1] P )) VP } 
= supj^A 1 ^ - if(F,d,p)A p V^*\ 

A>0 L J 

= (Hd/(p*K(H,d,p))) Hd/{p *- Hd \l - Hd/p*)(6/p) p *^ p *- Hd l 
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Consequently, 

liminfa-P'/^^logEexpj^a 217 ^^ [Lf(B H )] p dx^j *\ (2.11) 
> p- 1 {Hd/(p*K{H,d,p))) Hd ' { P*- Hd \l - Hd/p*)(e/p)P*/ < -P*~ Hd \ 

If this can be strengthened into equality with limits, then by Gartner-Ellis theorem, for any A > 0, 

Jim a-v*/( Hd ^ logpj (J [L X 1 (B H )] P dx^J > Xa^A 

f 1 Hd p* 

= - sup <X9- p (Hd/(p*K(H, d,p))) v^Hd (i-Hd/p*) {9 /p) \ 
= -p~ 1 K(H,d,p)\ I >''/( H V 

In particular, 

lim a- p * /(/Mp) logp( f \Lf(B H )V dx > a] = -p' 1 K{H,d,p). (2.12) 

The conjecture (|2.12p is partially supported by a recent result of Hu, Nualart and Song (Theorem 
1, |27| ) which states that when Hd < 1 and p = 2 



[Lf(B H )] 2 dx^ <C n {n\ 



nr^Hd n = 1)2 ,... 



for some C > 0. Indeed, a standard application of Chebyshev inequality and Stirling formula leads 
to the upper bound of the form 

lim sup a- 1/{Hd) logpj I [Lf(B H )] 2 dx > XaX < -I, 

where I is a positive constant. This rate of decay of tail probabilities is sharp by comparing it with 
(I2TTD for p = 2. 

In the case Hd > 1, /?^([0,t]<) can not be properly defined. On the other hand, this problem 
can be fixed in some cases by renormalization. For simplicity we consider the case p = 2. Hu and 
Nualart prove (Theorem 1, [26j) that for 1 < Hd < 3/2, the renormalized self-intersection local 
time formally given as 



7 H ([0,t] 2 <) = // S {B H (r) - B H (s)) drds 



'{0<r<s<t} 

-E ft 5 (B H (r)-B H (s))drds 

J J{0<r<s<t} 

exists with the scaling property 

i H {ml)=t 2 - Hd i H {m 2 <) (213) 

We also point that an earlier work by Rosen (|42j) in the special case d = 2. 



S 



Based on a similar but more heuristic reasoning, it seems plausible to expect that 

lim a- 1 /( ffd )logP( 7 /f ([0 ) lE) >a\ = -l^'^KlH, d, 2) (2.14) 

a— s>oo L J 

We refer the interested reader to Theorem 4, |27j for some exponential integrabilities established 
by Hu, Nualart and Song based on Clark-Ocone's formula. We leave these problems to the future 
investigation. 

Our large deviations estimates can be applied to obtain the law of the iterated logarithm. 



Theorem 2.5 When Hd < I, 

limsupt-^ Hd \loglogt)- Hd L° t {B H ) = 9(H,d)- Hd a.s. (2.15) 

t— >oo 



\\^n V t-^- Hd \\og\ogt)- Hd L° t {W H ) = e{H,dy Hd a.s. (2.16) 

t— >OD 

When Hd < p* , 

limsupt-P^-^^Ooglog^-^-^a^dO,^) = K(H,d,p)- Hd{p - 1) a.s. (2.17) 

t— >oo 



limsupt- p{1 - Hd/p *\loglogt)- Hd{p - 1) a H ([0,t} p ) = K(H,d,p)- Hd{p ~ 1) a.s. (2.18) 

t— >oo 

Even with the large deviations stated in Theorem 12.11 — [2^1 the proof of Theorem 12.51 appears to 
be highly non-trivial due to long-range dependency of the model. Here we mention some previous 
results given in Baraka and Mountford (|3f); Baraka, Mountford and Xiao ([!])• Using the large 
deviation estimate similar to (12. ip . Baraka, Mountford and Xiao were able to establish some laws 
of the iterated logarithm which describe the short term behaviors (as t — > + ) of the local times of 
fractional Brownian motions. As pointed out by Baraka and Mountford (p. 163, [3]), their method 
does not lead to the laws of the iterated logarithm of large time given in Theorem 12.51 

Theorem 12.51 will be proved in section [6j The proof of the lower bound appears to be highly non- 
trivial due to long-range dependency of the model. The approach relies on a quantified use of 
Cameron-Martin formula. 

Since all main theorems stated in this section have been known in the classic case H = 1/2 (see, 
e.g., [10J and [13J), we assume H ^ 1/2 in the remaining of the paper. 



3 Basic Tools 



In this section we provide some basic results that will be used in our proofs. We state them 
separately for a convenient reference. 
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3.1 Comparison of local times 



We will give general comparison results for local times for Gaussian processes. They are based 
on the standard Fourier analytic approach but go far beyond, motivated mainly by similar small 
deviation estimates. We start with an outline of the analytic method typically used in the study of 
local times for Gaussian processes, in particular on its the moments, see Berman [8] and Xiao |47| . 

For a fixed sample function and fixed time t > 0, the Fourier transform on space variable x £ M. d is 
the function of A £ R d , 

e iX - x L{t,x)dx= [ e iX - x ^ds. 



o 



Thus the local time L(t, x) can be expressed as the inverse Fourier transform: 

L(t,x) = -^— d [ e- iA - fe iXx ^dsdX. 

The m-th power of L(t, x) is 

L (*' x ) m = 7T"W I e~ ix -^ Afc f exp (i V X k ■ X(s k )) ds x ■ ■ ■ ds m dXi ■ ■ ■ d\ m . 

[ZlT) J R md J\0,t\ m 



Take the expected value under the sign of integration: the second exponential in the above integral 
is replaced by the joint characteristic function of X(si), • • • , X(s m ). In the Gaussian case, we obtain 

EL(t,x) m 

= to w / e" ix ' E ™=i Afc / exp f - -Var( V \ k ■ X(s k ))) da x ■ ■ ■ (2s m (fAi • • • dX m . 
K Z7r ) JR md J[0,t] m v 1 J 

Interchanging integration and applying the characteristic function inversion formula, we can get 
more explicit (but somewhat less useful) expression in terms of integration associated with det (EX (si)X(sj))' 
Estimates of the moments of local time L(t, x) thus depend on the rate of decrease to of det(EX (si)X (sj)) 
as Sj Sj — i — > for some j. Here in our approach, we have to make proper adjustment by approx- 
imating L(t, x). 

Consider now a random fields X(t) taking values in R , where t = (t±, . . . , t p ) £ (M. + ) p . For a fixed 
Borel set A C (M + ) p , recall that the local time formally given as 

L x (A,x)= [ 6 x {X(s))ds (3.1) 

is defined as the density of the occupation measure 

H A (B)= [ l B (X(s))ds BcM. d 

if ha(') is absolutely continuous with respect to the Lebesgue measure on 
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Given a non-degenerate Gaussian probability density h(x) on M. d and e > 0, the function h e (x) 
e~ d / 2 h(e~ 1 / 2 x) is also a probability density. Define the smoothed local time 



L x (A,x,e) = / h e (X(s) - x) ds. (3.2) 

J A 



Our first proposition provides moment comparison (13. 6p which can be viewed as analogy of Ander- 
son's inequality in the small ball analog: For independent Gaussian vectors X, Y, X symmetric, 

P(||X + y|| < e) <P(||X|| < e). 

See Li and Shao [M] for various application of this useful inequality. 

Proposition 3.1 Let A C (M + ) p be a fixed bounded Borel set. Let X(t) (t = (i x , . . . ,t p ) G (M + ) p ) 
be a zero-mean M. d -valued Gaussian random field with the local time L x (A,x) continuous in x £ W 1 . 
Assume that for every m = 1,2,... 

/ dsi • • • ds m [ d\i ■ ■ ■ d\ m exp { - -Var ( Y] X k ■ X(s k )) 1 < oo . (3.3) 
Then L x (A,0) € C m (i.e., finite m-th moment), with 



1 



EL X (A, 0) m = — — d I dsi ■ ■ ■ ds m I dAi • • ■ dA m (3.4) 



f 1 m 1 

exp i - -Var ( ^ A fc • X(s fc )) I 

^ fc=i J 



and 



lim E\L x (A,0,e) - L x (A,0)\ m = 0. (3.5) 

IfY(t) (t = (t\, . . . ,t p ) G (M + ) p y ) zs another zero-mean M. d -valued Gaussian random field indepen- 
dent of X(t) such that the local time L x+ y(A,x) of X(t) + Y(t) is continuous in x, then 

E [L X+Y (A, OH < E 0) m ] . (3.6) 

Proof: By Fourier inversion, we have from (13.21) 

MA, 0. „) - ^ I dAexp { - 1 (A • FA)} / 
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where V is the covariance matrix of Gaussian density h{x). Using Fubini theorem, 

1 



EL x (A,0,e) m = — — d I dsf-dam I dXi---dX r , 



exp I ~ Xk ' r ^ fe f Eex P \ ~i^2^k- X(s k ) \ 
^ k=l ' k=l ' 



(3.7) 



(27T) 



md 



dsi ■ ■ ■ ds m / dXi ■ ■ ■ dX Ti 



( m 1 f 1 m 1 

exp I ~ \ ^2 Xk ' TXk r exp i ~ 2 Var {^2 Xk ' x ( Sk ^) r • 



By monotonic convergence theorem, the right hand side converges to the right hand side of (13. 4h as 
e — > + . In particular, the family 

EL x (A0,e) m (e>0) 

is bounded for m = 1, 2, • • • . Consequently, this family is uniformly integrable for m = 1, 2, • • • . 
Therefore, (|3.4p and ([3.5p follow from the fact that Lx(A,0, e) converges to Lx(A, 0), which is led 
by the continuity of Lx(A, x). 



Finally, (I3.6P follows from the comparison 



/ tfei ■ • ■ ds m / tZAi • • • dA m exp I - ^Var ( V] A • (X(s k ) + Y(s k )) 1 
< / tfei • • ■ ds m / dAi • • • d\ m exp { - - Var (j^X- X(s k )) \. 

□ 

In certain situations we can also reverse bound in (|3.6p as a result of the Cameron-Martin Formula. 
In small ball setting, this is motivated by the Chen-Li's inequality [12J which can be used to estimate 
small ball probabilities under any norm via a relatively easier L2-norm estimate. See also the survey 
of Li and Shao |34| . Let X and Y be any two centered independent Gaussian random vectors in 
a separable Banach space B with norm || • ||. We use | • |^(x) to denote the inner product norm 
induced on by fj, = C(X). Then for any A > and e > 0, 

F(\\X + Y\\ < e) > F(\\X\\ < e) • Eexp^^YlSW' 

and 

^11 < e) > P(||X|| < Ae) • Eexp{-2- 1 A 2 |y|2 }. 



'/*(*)- 

Next we provide the local time counterpart of this inequality, which is crucial in our estimates. 
Suppose that the process X(t), t G [0, T], where T = (Ti, . . . ,T p ) G can be viewed as a 

Gaussian random vector in a separable Banach space B such that the evaluations x i— > x(t) are 
measurable (say B = C([0, T]; for concreteness) . Let %{X) denote the reproducing kernel 
Hilbert space (RKHS) of -X"(t), t G [0, T] equipped with the norm || • ||. Now we will make a crucial 
assumption that the independent process Y{t), t G [0, T] has almost all paths in %{X). 
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Proposition 3.2 In the above setting, under the assumptions of Proposition COl we have 

E [L X+Y (A, 0) m ] > E e -ill y H 2 E [L X (A, 0) m ] , (3.8) 
for every A C [0, T] and m £ N. 

Proof: Applying Lemma [3~W ii) . for g(x) = IJfcLi ^e(^(sfc)), x £ B, we get 

E[L x+Y (A,0,e) m } = / ds 1 ...ds m ETT/ le (X(s fc ) + y(s fc )) 
^ m fc =i 

„ m 

>Ee-^ n ' 2 / dsi E]T/i £ (X(s fc ))) = Ee^ll y ll 2 E [L X (A 0, e) m ] • 

jAm k=i 

Applying f|3. 5j) for both processes, X and X + Y, we get (|3.8p . □ 

3.2 RKHS of W H (t) and the remainder 

Let H £ (0, 1/2) U (1/2, 1) and recall decomposition (fL6)) : 

c^B H (t) = W H {t) + Z H (t), t>0, 
where the remainder process Z H {t) can be written as 

Z H (t)= {(t + sf-^-s^/^dBis), (3.9) 
•/ o 

with B(s) := B(-s), s > 0. Clearly, is a self-similar process with index H and the processes 

W (i) and Z H (t) are independent. In this section we develop a technique allowing us to treat 
sample paths of Z H (t) as, essentially, elements of the reproducing kernel Hilbert space (RKHS) of 
W H (t). 

The RKHS H[0,T] of the the Riemann-Liouville process {W H (t)} yi with index H > 0, viewed 
as a random element in C[0,T], follows standard theory of RKHS, see [33 J and [6]. Van der Vaart 
and van Zanten [44 [ Lemma 10.2] proved that 

H[0,T] = /£ +1/2 (L 2 [0,T]), (3.10) 

where 

Io + f(t) = f ^ ) j\t-sT- 1 f(s)ds, te[0,T] (3.11) 
is the Riemann-Liouville fractional integral of order a > 0; for a = 0, := /• 

Proposition 3.3 {Z H (t)} t> has C 00 -sample paths a.s. for any a > 0. However, for every T > 

P({^)} te[0 , T] eH[0,T])=0. 
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Proof: Formal ra-tuple differentiation of Z H (t) gives 

|V(t) = f[(H - 2 -^l) jH(t + s) B -<**W dB(s), t > 0. 



k=l 

The right hand side is a well-defined Gaussian process with locally square integrable sample paths. 
By consecutive integration of this process over [a,t] we prove that {Z H (t)} f>a has C^™ -1 ) -sample 
paths, n > 1, which proves the first part of the proposition. 

To prove the second part observe that 

Z H (t) = (I^ 1/2 V H )(t), t > 0, (3.12) 

where V H (t) is a Gaussian process given by 

H - 1/2 [°° t- H -V*u H -V 2 



Direct computation gives E (V^(i)) = Ci 1 , where C depends only on H. Hence E|| H^rp T j = 

oo but E|| V r ^||i 1 [o ) r] < oo. Combining the fact that I^ 1 ^ 2 is one-to-one on Li[0,T] (see |43[ 
Theorem 2.4]) with (pTIgJ) and (l3~TUI) we get 

p W] eH[0,n) =P({^W} te[0 , r] eZa[0,r|) =0, 
where the last equality follows from a zero-one law and integrability of Gaussian noms. □ 

Direct verification whether a given function belongs to H[0, T] can be difficult. Therefore, we give 
below a simple to check sufficient condition. Let j4C™[0, T] denote the space of functions / which 
have continuous derivatives up to order m— 1 on [0,T], with /( m_1 ) absolutely continuous on [0,T], 
and /( m ) G L 2 [0,T], m G N. 

Proposition 3.4 Let m = \H + 1/2]. If f G ^C^[0,T] is suc/i that / (fe) (0) = for < k < m, 

then f G H[0,T] and 

ll/llH[o,n = k H \K+ iH+1/2) f^h^Tl , (3.13) 

u?/iere fc H = T(H + 1/2) -1 . 

Proof: By our assumption / = /^/ (m) , where /( m ) G L 2 [0,T]. Put g = I ™-( H + 1 / 2 ) f(m) _ gince 
the operators of fractional integration {Iq + : a > 0} form a strongly continuous semigroup on 
L 2 [0, 1] (see 03 Theorem 2.6]), we get that G L 2 [0,T] and 

j-H+l/2 _ T H+l/2 ( T m-(H +1/2) ,( m )\ _ T m Am) _ f 

1 o+ 9 — 1 o+ y 1 o+ J ) — J o+J — J ■ 

In view of (f3TTU]) . / G M T and from gH Lemma 10.2] 



7 ii ii , ii T m-(H+l/2) Arn)\\ 

H[0,T] = KH\\g\\L 2 [0,T] = k H\H 0+ J '\\l 2 [0,T] 



□ 



The remainder Z H is not in H[0, T] by Proposition [3731 The next result shows the way to circumvent 
this problem, which is crucial to our technique. 
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Proposition 3.5 For any a > there is a Gaussian process |Z^(t)} t>Q such that 

(i) Z**(t) = Z H {t) for all t > a; 
(ii) for any T > 



Proof: First consider H G (0, A), so that m = \H + 1/2] = 1. Define 



/At, 

\z H (t), 



0<t<a 
t > a. 



where A = a 1 Z H (a). Since Z^(t) has paths in AC^lO, T] (see the first part of Proposition 13.3 
and Z^ (0) = 0, (ii) holds by Proposition 13.41 

Now we consider € (|, 1), so that m = \H + 1/2] = 2. Define 



Bit 2 + B 2 t 3 , 0<t<a 
Z H (t), t > a 







where Bi = 3a" 2 Z H (a) - a" 1 Z H (a), 5 2 = -2ar z Z H {a) + a~ 2 Z H (a), and Z H (t) := Q^Z H (t). As 
in the previous case, part (ii) follows by Proposition 13.41 Indeed, Z^(t) has paths in ^IClfOjT], 
Z* (0) = 0, and (0) = 0. □ 

The above method of modifying of Z H in a neighborhood of will also be used in Section [6] for 
other processes and the H[0, T]-norm of such modifications will to be estimated. For this purpose 
the next lemma will be useful. 



Lemma 3.6 Let m = \H + 1/2]. Iff G AC^[0,T] and / (fc) (0) = for < k < m, then for every 
a G (0,T) 



2 

H[0,T] 



< c f (T2m - 2 H_ a2m ^ )||/M||L[M+ /" T / 
I J a J a 



(t- s) m ~ H - 3/2 f( m \s)ds 



dt 



where C depends only on H . 



Proof: Put k = m - {H + 1/2). In view of (f37T3|) we get 

Ih[0,T] = ^ll J 0+(/ (m)l [0,a] + f im)l [a,T])\\L 2 [0,T] 

< 24ll'oV[O l0 ]|li 2[ 0,T] ll/ (m) HL[0,a] +2k 2 H \\I«Af (m) ha,T])\\l 



[0,T] 



< C(T : 



2m-2H 2m-2H \u f (m) \\2 

)\\J llLoo[0,a] 



a 



+ 2k 2 H 



T 



T 



(t- s r- i f ( - m \ s )ds 



dt. 



□ 
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3.3 Technical lemmas 



The following auxiliary results and formulas are used in the proofs of main theorems. They are 
given here for a convenient reference. 

Lemma 3.7 Let fi be a centered Gaussian measure in a separable Banach space B. Let g : B i— )■ M + 
be a measurable function. Then 

(i) if {x E B : g(x) > t} is symmetric and convex for every t > 0, then for every y E B 

g(x + y) fi(dx) < / g{x)n(dx); 

B J B 

(ii) if g is symmetric (g(—x) = g{x), x E B), then for every y in the RKHS rl^ of /i 

J^g{x + y) n{dx) > exp |-^||y||^| g(x) n(dx), 
where ||y|L denotes the norm in 

Proof: Part (i) follows from Anderson's inequality 

p poo 

/ g(x + y) fi(dx) = / fi{x E B : g(x + y) > t} dt 
Jb Jo 

poo p 

< / fi{x E B : g{x) >t}dt= / g(x) /J,(dx). 
Jo Jb 



Part (ii) uses Cameron-Martin formula and the convexity of exponential function 

'l(-r) exp <J {x.y),, - 

is Jb 



g(x + y) n(dx) = I g{x) exp { (x, y)^ - ~ \\y ||£ \ fi(dx) 



: ^^5 f (^)exp|(x,y) Ai -i| 



y\\l )■ n(dx) 



1 r , , r i.. , 2 



+ 2 J B 9(x)exp <j -{x,y)n - -\\yWp, ) n(dx) 



>exp|-^||y||^| J^g(x)fi(dx) 



□ 



The next lemma is well-known and goes back at least to 1950s in equivalent forms, see Anderson 
[TJ p. 42], Berman (7J p. 293], and [8j p. 71]. The basic fact is that conditional distribution 
of Xk given all the Xi, 1 < i < k is a univariate Gaussian distribution with (conditional) mean 
E(Afc|Ai, . . . ,Xf t _i) and (conditional) variance 

det(Cov(A 1 , . . . , X k ))/ det(Cov(X 1 , . . . , X k _{j) 

for 1 < k < m. 
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Lemma 3.8 Let (X\, . . . ,X m ) be a mean-zero Gaussian random vector. Then 

det(Cov(Xi, . . . , X m )) = Var(Xi)Var(X 2 | Xt) ■ ■ ■ Var(X m | X m -i, ■ ■ ■ , X x ). 

Let B H (t) be given by its moving average representation ([1.4p . By the deconvolution formula of 
Pipiras and Taqqu [40J we also have 



B{t) = c* H ! Ut - s) 

J — oo ^ 



1/2-H 



-s)f~ H ) dB«(s), 



(3.14) 



where c* H = {chF(H + l/2)T(3/2 — H)} 1 and the integral with respect to B H (t) is well-defined 
in the L 2 -sense. It follows from (jl.4p and (|3. 14j) that for every t £ R 



H-i 



T t ■= a {B H (s); -oo < s < t] = a{B(s); -oo < s < t}, 
where the second equality holds modulo sets of probability zero. Then for every s < t 

E(B H (t)\T s ) = c H f ((* 
If d = 1, then for every s < t 

Var (B H (t) | T s ) = e{ \b h {t) - E(B H (t) \ T s ) ~ I 



dB(u) 



(3.15) 



(3.16) 



J\t-u) H -*dB(u)\r a \ 

c 2 h j\t-u? H - 1 du = ^L{t-s\ 



2H 



(3.17) 



For the reader's convenience we also quote the following lemma due to Konig and Morters, 
Lemma2.3]. 



Lemma 3.9 Let Y > be a random variable and let 7 > 0. If 



lim — log 

m— >oo rn 



!Y7 



-EY r ' 



for some k£R, then 



lim ~tk 1o § p { y >y} = -7 e " 



-«/7 



(3.18) 
(3.19) 



4 Large deviations for local times 

4.1 Proof of Theorem 12.11 — superadditivity argument 

In the light of Lemma 13.91 it is enough to show that the limit in f|3. 18j) exists for Y = L\{B H ) 
and for 7 = Hd. We will prove it by a superadditivity argument. Let r be an exponential time 
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independent of B H (t). We will first show that for any integer m, n > 1, 



E 



/ r>H\m+n 



>\ m + n )E 

m 



L° T (B H ) m E L U T (B H ) 



(4.1) 



Let t > be fixed. Notice that by Theorem 17. 1[ the Gaussian process B (t) satisfies the condition 
(|3.3p posted in Proposition 13.11 By (|3.4p . therefore, 



E 



L° t {B H ) m 



,„ \ mH [ dsf-dsm [ dXi ■ ■ ■ a?A m exp \ - ^Var fV] \ k ■ B H (s k )) \ 
(2ir) md J [0)t]m J(Rd)m { 2 V J J 



io. 



dsi • • • ds r , 



(2tt)^ y [0it] 

where Bq (t) is 1-dimensional fractional Brownian motion. 
By integration with respect to Gaussian measures 



f ( 1 m 1 

/ dXt ■ ■ ■ d\ m exp < - ^Var X k B^(s k )) I 



/ rfAi • • • dA m exp | - ^Var ( X k B^(s k )) j 
= (2vr) m / 2 det {Cov (b$ ( Sl ), • • • , 5^ ( Sm )) }~ V2 . 



Therefore, 



E 



0/ r>H\m 



L» t (B 



ml 



ds! • • • ds m det {Cov (b^(si), • • • , Bg(s m f) } 



-d/2 



( 27r )md/2 y [()>tK 

In (I4.2j) and elsewhere, for any ^4 C M + and an integer to > 1, we define 

A™ = {(si,--- ,s w ) E si<---<s m }. 



Put 



and «4(si, 



,Sk) = v{Bo(si), ■ ■ ■ ,Bo(s k )}, fe = 1, • • • ,to, 
, Sfc) = {0, f2} when = 0. By Lemma |3.8[ 



E 



( r>H\m 



dsi ■ ■ ■ ds m cp m (si, ■ ■ ■ ,s m ), 



where 



( 27r)m d/2 
m 

<Pm(si, ■ ■ ■ , %) = II Var ( 5 o H (sfe)l^(si), ■ • • , B?(s fc . 

k=l 



-d/2 



(4.2) 
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with the convention that the first term is Var {Bq (si)) for k = 1. We are ready to establish (|4.ip . 
Let m, n > 1 be integers. Then, for any si < ■ ■ ■ < s n + m and n + 1 < k < n + to, 



Var {Bg(s k )\Bg( 8l ), • • • , 5 % fc _i 

= Var (<(s fc ) - • • • , 

= Var (B(?(8 k ) - B*(s n )\BX( Sl ), • • • , B* \s n ), B* (s n+1 ) - B*(s n ), 

■■■ ,<( Sfc _!)-<( Sn ) 

< Var to(s fc ) - B*(s n )\B*(s n+1 ) - B?(s n ), ■ • • , - B^(s n 



= Var (#^(s fe - s n )|B^(s n+ i - s n ), • • • ,Bff(s k -i - s n ) 
where the last step follows from the stationarity of increments. Thus 

) > <Pn(si, • • ■ , S n )ip m (s n+1 - 



Notice that from (IP 



E 



/ TjH\m 



ml 



E 



dsi • • • ds m ip m (si, ■ 



ml 



(2ir) md / 2 

m\ 
(2ir) md / 2 



E 



1 



Sm<T 



dsi - ■ ■ ds m ip m (si, - ■ ■ ,s m ) 



Sl<-<Sm 

dsi ■ ■ ■ ds m tp m (si, ■•■ , 



Consequently, 



E 



( r>H\n+m 



L»(B») 



> 



(n + m)\ 

' /TTp _|_ \ Ti -|- 77T 

fn + m)! 



dsi---ds n+m (f n (si, 



Sn e 



( 27r )(n+m)d/2 y (R+) „ 

X ^m(Sn+l j ^n+m Sn)& n + m 

(n + m)\ f , . _ , 

dsi---ds n v? n (si,--- ,s n )e n 



dti • • • di m <Pm(ti, • • • , t m )e 



( 27r )(n+m)d/2 _/ (R+) „ 



n + m 



■/?? 



E 



L»{B 



E 



L° T (B H ) m 



We proved relation (|4.ip that says that the sequence m \— > log 
By Fekete's lemma the limit 

lim — log | -^E 



rE 



ml 



(4.3) 



L®(B H ) m I is super-additive. 



TO! 



Ofr>H\m 



L»(B 



1 / 1 
sup — log — -E 



f r>H\m 



I? T {B 



(4.4) 
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exists, possibly as an extended number. By the scaling property (jl.ip . 



E 



L° T (B H ) m 



E 



_(1-Hd)m 



E 



Li(B H ) m = T (1 + (1 — Hd)m) E V{{B H ) 



0[ r>H\m 



From (|4.4p and Stirling's formula we get 



lim — log 



7 E 



(ml) Hd 

Applying Lemma 13.91 we establish (|2.ip with 

6(H, d) = Hd(l - Hd)- 1+1 / Hd L~ 1 / Hd . 



To obtain (|2,2p and complete the proof it is enough to show that 



(2 7 r)- a!/2 r(l -Hd)<L< (H-^cjj) T(l - Hd) . 



-d/2 



By gTj 



-E 



ml 



0( r>H\m 



Ut(B 



> {EL° T (B H )} m = {(2vr)- d / 2 r(l - Hd)}™ , 
where the equality comes from (|4.3p (for m = 1). This proves the lower bound in (14. 7\i . 
To prove the upper bound, we first notice that 

Var (Bg (s k ) | B^( Sl ), ■ ■ ■ (s^)) > Var (b£ (s k )\ B (s),s < s k ^ 

r 2 

— H ( a _ \2H 

~2H^ Sk ~ Sk ~ l > ' 
where we used (|3.17p . Hence the function (p defined above satisfies, with sq = 0, 



Vm(si)--- ,s m )< (2H / c 2 H ) md/2 JJ(sfc - Sfc-i 



-Hd 



k=l 



and by (|43|1 . 



md/2 



E 



L°JB H ) m 



<ml dsi ■ ■ ■ ds m TT(sfe - Sfc-i) Hd e" 



fc=l 
m 



r poo i m 

ill J t- Hd e- t dt\ =m\T(l-Hd) r 

This establishes (|4.7p and completes the proof. □ 

4.2 Proof of Theorem 12.21 — comparison argument 

First we note that 

L° t {c- H X B H ) = c H L° (B») . 
Thus, from the decomposition (|1.6p and f|3 .6[) for every m £ N, 



< E 
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To prove a reverse inequality (up to a multiplicative constant) we use notation (|3,ip . Fix a £ (0, 1) 
and let let Z^(t), t > be the process specified in Proposition 13.51 that is also independent of 
W H (t), t>0.We have 

c d H L\ (B H ) =L c -x BB ([0,l],0) >L c _i BH ([ a ,l],0) =L wH+Z H([a,l},0). 

Thus, by fj3.8|) we get 



cg d E Ll(B H ) n 



>E 



L wH+Z H([a,l],0) m \ > K a E[L wH ([a,l],0) m ] 
=K a E[(Ll(W H )-L° a (W H )) m ' 
>K a {e [L\(W H ) m ] 1/m - E [L° a (W H ) m ] 1/m } m 



K a (l-a 



l-Hd 



E [L$(W 

where the last equality uses self-similarity (11 .3^ and K a = Eexp { — ^||Z^|| 2 } . Here ||^<^|| < oo a.s. 
is the RKHS norm associated with {W \t)}te[o,i] an d computed for paths of {^^(i)}te[o,ll ■ This 
together with (14, IIP yields 



c^E L\{B H ) m <E L\{W H ) 



-0 (urH\™ 



< K7 1 ( 1 - a 



A-Hd 



c^E L?(B H ) m 



Applying the limit as in (|4.5p to both sides and then passing a — >• gives 



lim — log 



1 



>!) 



\\Hd 



E 



L?(l^ H ) m ] ^ = log {4(1 - Hd)-( 1 - Hd >Ly 



Therefore, by Lemma 13.91 the limit in fj2.5j) exists and 0(H,d) = c^^ H 9(H,d) by (|4.6 



□ 



5 Large deviations for intersection local times 



5.1 Proof of Theorem 12.31 — subadditivity argument 

Let (A) be defined analogously to (jl.9p by 



af(A)= ! [ f\ Pe (Wf(sj)-x)ds 1 ---ds p dx, 

JR d J A - =l 

where p e is as in (jl.lUp . We will first prove the subadditivity property: for every m, n £ N, 



E 



af ([0,n] x ••• x [0,7-p])' 



m+n 



(5.1) 



f m + n 

< 

V m 



E 



5?' 



([0,n] x ••• x [0,r p ]) m Maf ([0,n] x ••• x [0,r p ]) n 
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where t\ , . . . , T p are iid exponential random variables with mean 1 and independent of Wf(t), W^(t). 
Indeed, since 



f P m . t . 

df ([0,ii] x • • • x [0,t p ]) m = / dxi ■ ■ ■ dx m TT TT / Pe{Wf \s j)k ) - x k ) ds jjk , 

J(R d ) m 7L -, . i ./0 



3=1 fe=l ' 



we can write 
E 



df ([0,ti] x • • • x [0,T p ]) m+n = I dxi ■ ■ ■ dx m+n £(xi, . . . ,x m+n ) p , (5.2) 

J J(R d ) m + n 



where 

/>oo p m-\-n 

£(xi,... ,x m+n ) = dte~* dsi ■ ■ ■ ds m+n E TT p e {W H (s k ) - x k ). 

Jo J[o,t]™+* £J[ 

Let 

A = {(si, • • • ,s m+n ) G [0,i] m+n : max{si, . . .,s m } < min{s m+ i, . . .,s m+n }} . 

There are exactly ( m ^ n ) permutations Uj of {1, ... ,m + n} such that Ui (J i~ 1 A = [0, t] m+n and 
a^ 1 Dt are disjoint modulo sets of measure zero (here, a(si, . . . , s m + n ) := {s a m, ■ ■ ■ , s CT (m+n)))- 
Therefore, 



/ dsi • • • ds m+n E TT p e (W H (s k ) - x k ) 
J%t] m+n fc=i 

„ m+n 
» i lD ' fc=l 

dsi • • -ds m+n E p t (W H (s k ) - x ai ( k )), 
i Jt k=i 

which gives by Holder inequality 

£{xi, . . . ,x m+n ) p = ly^l dte~ l I dsi ■ ■ ■ ds m+n E |J_ p € (W n (s k ) - x a .^) 



< 



m + n 
rn 



(poo p m + n 

f = \ E / dte- f d Sl --- ds m+n E TT p e {W H (s k ) 

{i J ° Jd * fc=l 

\p-l ( roo r m + n Y 

) J2\j Q dt e_< j D dsi • • • dsm + n e n ^ - x -m) i 



Substituting into fj5.2|) yields 

p-l 



E 



f([0,r 1 ]x...x[0,r p ]) m+ "l < + W (fai-cfa^n 

{/■oo /■ m + n 

y ^ e "*y dsi---ds m+n E JJ p £ (W H (s fc ) -x CT . (fc) ) 



) dx! ■ ■ ■ dx m+r 



^ dt e-* y^ dsi • • • ds m+n E JJ p £ (W H (s fe ) - x k ) | . 
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Since the last integrand can be written as 

\ p 



( poo p m + n Y p 

I die' 1 dsi---ds m+n M TT Pe (W H {s k ) -x k ) \ = 

[JO JD t ^ J J(R+)P 

p m+n p 

x / ds^i ■ ■ ■ ds jtm+ n) E Yl Y[Pe{Wf(8 jth ) - x k ), 



after integrating with respect to an, ... , x m + n we get 



E 



of ([0,n] x ••• x [0,r p ]) m+n < ( m + n ) P / p dt 1 ...dt p e-^+^ (5.3) 



m 

m+n 



where 



/ ( n dS ^ ' ' ' ds hm+n) E ]J QeiWfis^k), W* (a p , k )) , 

f P 

g e {yi,--- ,y P ) - = / Y\pe{yj-x)dx 



(5.4) 



(27re) -*-i)/2 p - rf /2 exp { _ 1 £ | y . _ y|2}) 



2e • i 



and y : = p 1 X]f=i ^« f° r 2/1 > • • • > Vp ^ ^ • Moreover 

p m+n 

Q dsj 5 i • • • ds j)7n+n ) E J| 

lD H X-xD tp fe=1 



/ ( J! ds jA ■ ■ ■ d Sj , m+n ) E H g € (Wf(s ltk ), W*{a Ptk )) (5.5) 

JD H X-xD tp V J=1 fe=1 
r p r p 

/ ( TT ds il " ' ds j,m) / (T\ds j>m+1 ---dsj ;m+n ) 

J[o,t] m K fj[ ' y[o,t-s*] n v J= i ' 

m m+n 

xEl[g £ (Wi I (s 1>k ),...,W p H (s p>k )) If ge[Wf(8t + 8 1)k ),-" ,W*(8* p + s P)k )), 



where 

t = (ti,--- ,t p ), s* = (si,-- - ,s* p ), 

and 

s* = max{sj fc : 1 < k < m}. 

Assuming that W^(t) are given by (|1.2p with independent Brownian motions Bj(t), define J-~ s * 
<j{i?j(nj) : tij < Sj, j = 1, . . . Put also 



YJ- ( s * , s) = / (s*+s- 5 dBj (u) and Z(s* , s) = / (s* + 

is* " ' io 



s - u) H -^dBj(u), 
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so that Wj^Sj + s) = Yj(s*, s) + Zj(sj, s). The last expectation can be written as 



x E 



fc=i 

m+n 

[ 5,(^(4,^) + zf (4, Sljfc ),.-. ,y/( s ;, Sp)fe ) + zf( s ; jSp)fc ))|^]} 

fc=m+l 

m m+n 

<^[\\g e (wF{s lM ),...,W^(s p , k ))y[ H ge{Y 1 H (8* li8ltk ) i --- ,Y p H (s* p , s Pjk )) 

k=l k=m+l 
m m+n 

= E[Hg e {W 1 H (s hk ),...,W p H (s p , k ))]E[ H g e {Wf(si, k ),--- ,W?(s Ptk )) 



k=l k=m+l 

where the inequality follows from Lemma l3,7f i) (see the evaluation of g e in (|5.4j) and the positive 
quadratic form associated with it) and the last equality follows from 

[Yi(s*,si }k ), . . . ,Y p (s*, s Pjk )) = (W 1 H (s hk ),...,W p H (s p , k )). 

Combining the above bound with (|5.5p and then with (|5.3p we obtain 



E 



a 



f([0 1 n]x...x[0 1 r p ]) m ^ < 



m+nl . (m + n 



in 



dh ■ ■ ■ dt p e-( tl+ - +t ^ 



V m 

/ ( n ds i+ • • • ds ^Y n 9e{w?(s lik ), wf(s P}k )) 

J[o,t] m j=1 k=1 

m+n 

ds jtm+n ^E H g e (Wf( Sltk ),--- ,W p H (s p , k )) 

k=m+l 

p m 

( H ds jA ■ ■ • c% m )E 11 g e (wf( Sl>k ), W p H (s pM )) 



I I rfSj,m+l ' ' ' 
[0,t- 8 *]» v ^ =1 



m + n 
m 



x e 



-(*!+•••+*;) 



> j=l k=l 

dt x ■ ■ ■ dt v e-^- s *^ + - + ^- s p^ 

m+n 



[s*,oo]p 



^[0,t— S*]™ -J 7-_~ I 1 



fc=m+l 



m + n 



m 



E 



[of ([0,ri] x ••• x [0,r p ]) m ] l[of ([0, n ] x ••• x [0,r p ]) n 



where in the last equality we use 

e -(s*+-+ s £) _ 



-{ti+-+tp) 



Y\ l[s*,t]( s l,fc; • • • > s p,k) dh ■ ■ ■ dt % 



-) P k=l 

and the definition of g t in (|5.4I) . The subadditivity ([5.ip is thus proved for any e > 0. 

Now we would like to take e — > + on the both sides of (|5.ip in an attempt to establish 

Ea H ([0,r 1 ]x...x[0,r p ]) m+n 
m + n x 



< 



Ea H ([0,n] x ••• x [0,r p ]) m Ed H ([0,r 1 ] x ••• x [0,r p ]) n . 



(5.6) 
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To this end we need to show that for any m > 1, a H ([0, t±] 
and 



x • • • x 



[0,r p ]) is indeed in C m (n,A,P) 



lim E 



a 



([0,n] x ••• x [0,T p ]) m = E a H ([0,n] x ••• x [0,r p ]) 



(5.7) 



Indeed, using ([5.1 J) repeatedly we have that 



Notice that 



E 



E 



af ([0,ti] x ••• x [0,r p ]) m < (m!fl of ([0,n] x ••• x [0,r p ]) 



([0,rx] x ••• x [0,r p ]) 



e _ *Ep e (W ff (t)-a:)(ft 



e * / p e (y - x)p t * (y) dy 



dx 



where i* = {2Hy l t 2H and the last step follows from the easy-to-check fact that W H {t) ~ JV(0, (2H)- l t 2H I d ). 
By Jensen inequality, the right hand side is less than or equal to 



Pe(y - x) 



Pt*(y)dy 



dydx 



roo poo r 

/ ■■•/ dt 1 ---dt p e~^ + - +t ^ / X\ Pt *{x)d. 

Jo Jo J^ d J=[ 3 



JO 

<Z(p-l)/2 f°° f°° 



H/ir 



o Jo 



-(h+-+t p ) 



(e n 4 H y d/2 dh-dt 



j=l l<k=£j<p 

where the last step follows from a routine Gaussian integration. 
By arithmetic-geometric mean inequality, 



It n *i h s rr n * 

F 3=1 l<k^j<p j=l l<k^j<p 



2H/p 



JJ f 2ff(p-l)/p 



So we have 



E 



of ([0,n] x ••• x [0,r p ]) < (H/tv 



d(p-l)/2 



P 



-rf/2 



t - Hd{jp -l)/ Pe -t dt 



d(p-l)/2 ... 

ff/Tr) p" d/2 r(l - Hd/p*f. 



Summarizing our computation, we obtain 
(m!)- p E[af ([0,ri 



(5.8) 



By Theorem 17. 1[ the process 



x H (ti, ••• = (wf (ti) - wf (t 2 ), • • • , w^iCtp-i) - w?(tp)) 
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satisfies the condition fj3.3j) with A = [0, t] = [0, t\] x ■ ■ ■ [0, t p ] for any t\, ■ ■ ■ , t p > and 
df ([0,ti] x ••• x [0,i p ]) 



/ ^(Wf^O-Wf^),-- - ,W p H _ 1 (s p - 1 )-W p H (s p ))ds 1 ---d Sp 
J\o,t] 



p — 1 p— 1 



where 

■ /Rd i=i fc=j 

is a non-degenerate normal density on M^ -1 ). By Proposition 13, 1\ a^([0,ti] 
C m {Q,A,F) and 



x • • • x 



lim E 



5f ([0,ii] x ••• x [0,i p ]) m ] =E[a H ([0,ti] x ••• x [0,t p ]) 



[o,t p ]) e 

(5.9) 



In addition, by the representation (|3.7p one can see that for any e' < e, 



E 



fif ([0,ti] x ••• x [0,g) m <E dff([0,ti] x ••• x [0,t p ]) 



Thus, (|5.7p follows from monotonic convergence theorem and the identities 



E 



a 



f([0 > ri]x...x[0,r p ])' 



e -(tiH — Wp)]£ 



([o,ti] x ••• x [o,g) ? 



cZ£i • • • dtr 



and 



E 



a H {[0,n] x ... x [0,r p ]) m l = f e~^ + - +t ^E\d H {[0,h} x ••• x [0,t p ]) 

J J(R+)P 1 



(5.10) 



l • • • 



dt p . 



Further, by (|5.8|) we obtain the bound 

(m!)- p E[d i/ ([0,r 1 ] x ••• x [0,r p ]) m ] < r^/7r) d(P_1)/ V d/2 r(l-Hd/p*) p 



(5.11) 



The inequality (|5.6p implies that the sequence m i— >■ log ((m!) p Ed^([0, Ti] x 
sub-additive. Hence the limit 

lim — log((m!r p Ed H ([0,ri] x ••• x [0,rJ) m ) =c(H,d,p) 
m->oo 777, 

exists, possibly as an extended number. Further, by (|5.1ip 

\d(p-l)/2 _ 



x [0,t p ]) ) is 



c(H,d,p) < log |(if/7r) d(P lV V d/2 r(l - ^d/p*) p | . 



(5.12) 



(5.13) 



Now we will deduce the moments behavior of a ([0, l] p ). 
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Notice that t* = min{ri, • • • , r p } is an exponential time with parameter p. 

Ea H ([0,n] x ••• x [0,T p ]) m >m H {[0,n]P) m 
= Eri P ~ Hd(p ~ 1))m Ea H ([0, lf) m 

= p-(p-Hd(p-i))m r ^ + ( p _ Hd(<p _ i)) m ) Efi tf ([o, l]P) m . 

By Stirling's formula, 

limsup — \ogUm\y Hd(p ^Ea H {[Q,l] p ) m \ < c(H, d,p) - p{l - Hd/p*) log(l - Hd/p*). 
On the other hand, for every t%, . . . , t p > 0, 

Eaf ([0,ti] x ••• x [0,g) m 

„ p r m 

= dxf-dx m T] / rfsi • • • ds m ETTp e (ty H (s fc ) - 



i=i- y [°'* J ] m fc=i 

P ( „ / . m \P\1/P 



TT W dxi-dx m ( [ ds l ---ds m E\\p e {W H {s k ) - x k ) ) 1 



< 



Vp 



Letting e — > + , from fj5.9|) we get 

Ea H {[0,h) x • • • x [0,t p ]) m < f[ {Ea ff ([0,t,r) m } 1/p = Ea H ([0, l] p ) m • f[t^ 1 - Hd/p ' ) , 

3=1 3=1 

where the last equality uses self-similarity (|1.12p . Hence 

Ea H ([0,n} x ••• x [0,T p ]) m (5.14) 
= / dh ■ ■ ■ dt p e~ {tl+ - +tp) m H ([Q,ti] x ••• x [0,t p }) rn 

<Ed m ([0,lf) m I dt 1 ---dt p e-^ + - +t ^(h---t p ) m ^- Hd ^ 

J(R+)p 



Ea H ([0, l} p ) m r( 1 + m(l - Hd/p*)) P . 



By Stirling's formula again, 

liminf — log! (m\)- Hd ^Ea H ([0,l] p ) m \ > c(H, d,p) - p(l - Hd/p*) log(l - Hd/p*). 



m 

We have shown that 



lim — log\(m!)- Hd (P-VEa H ([0,l] p ) m \ = C(H,d,p), (5.15) 
%— >oc ml J 

C(H,d,p) = c(H,d,p) -p(l- Hd/p*) log(l - #d/p*) (5.16) 



m— i>oo 771 

where by f|5. 13[) . 



< log | {h/^ lV V d/2 r(l - Hd/p*) p (l - Hd/p*)- p( - 1 - Hd/p * ) ^ 
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On the other hand, let a H (A) be the intersection local time generated by c H Bi(t), • • • , CjfBjf(t). 
We have that 

a H (A) = 4 p - 1] a H (A), A c (R + f. 
In view of the decomposition (jl.6p . by Proposition 13 . 1 1 we have that 



E 



([o,if) 



> E 



a H (\0A]P) 



d{p-l)m w 



a* ([0,1]") 



It follows from (|5.24p below that 

d 

c( H , d, P ) > pio g |4 /p * (i - Hd/ P *r {1 - Hd/ ^ 2p * (1 + t 

Applying Lemma 13.91 leads the first conclusion (|2.7p of our theorem with 

C(H,d,p) 



K(H,d,p) = Hd(p-l)exp 



Hd(p 



2L\ 



and therefore the bounds given in ([2.8p follows from (j5. 16[) and (|5.19p . 



(5.17) 
(5.18) 



(5.19) 



□ 



5.2 Proof of Theorem 12.41 — comparison argument 

In connection to (|5.15p . we first show that 

lim — log((m\)- Hd( - p ^Ka H ([0,l} p ) m \ = C(H,d,p) - d(p - l)logc#, (5.20) 

The upper bound follows immediately from (|5.15p and the comparison ([5.18p . To establish the 
lower bound, we once again consider the intersection local time a H (A) generated by the normalized 
fractional Brownian motions 

B?(t) = c- H l B?(tl • • • , B?(t) = c- H l B*{t). 

For any e > 0, define 

af(A)= f [ T\p e (Bf(s j )-x)ds 1 ---dspdx, (5.21) 
Let < 5 < 1 be a small but fixed number. Notice 

„ m 

Eof ([0,lf) m >Eaf([(5,lF) m = / ds 1 ---ds m ET[gjB^( Sljk ),--- ,B^(s Ptk )) 

where g e (%i, • • • , x p ) is defined by f|5.4[) and we adopt the notation Sk = (si,fc, ■ ■ ■ s p ,fc)- 

Consider (Wjf^ii), ■ ■ ■ , W p (tpf) (t = (ti,--- ,t p ) £ [0, l] p ) as a Gaussian random variable tak- 
ing values in the Banach space (8>? =1 C{[0, l] p ,lR rf }. Then the reproducing kernel Hilbert space of 

(Wf(ti), • • • , W*{tpj) is i% = ® P j=1 H W - For each (/^x) • • • , f p (t)) G F w 

I|(/i(*i)--- ) /pW)HL = Ei^H^ 

i=i 
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where 1 1 • 1 1 h w is the reproducing kernel Hilbert norm of H\\r . 

Let Z^(t), • • • , Zj^(i) be the processes constructed in Lemma[33](with a = 5) by Z^ (i), • • • , (t), 
respectively. For each (si, • • • , s m ) S [5, l]P) m by the decomposition (|1.6p we have 



k=l 

m 

= E [] 5e(^( Sljfc ) + Zg^fc), • • • , Wf (s Pifc ) + Zf p (s M )) 



fc=l 



Fixed (si, • • • , s m ) 6 [S, lf) m . Applying LemmaE^ii) to the functional g(f u ■■■ , f p ) on <8? =1 C{ [0, 
defined by 

-,/ P ) = I]ft (/iKfe), • • • , /p(« P ,fc)) (/i, ■■-,/,»)€ ®? =1 C{ [0, If, R d }, 
fc=i 

then the right hand side is greater than 

(Eexp { - \\\Zf • • • , <) 

i \ p m 

:exp{--||Zf||^}j En^(wf(ai,*),'" >W? (***))■ 



Summarizing our estimate, we have 
Eaf ([0,lf) m 

Eexp{--||Zf||^}) / <fai---t^En^(W 1 H (ai,fc),--- 3 W^(s Plfc ) 

L 2 V ^([5,i]f) m fe =1 v 

= (Eexp{-i||Zf||^}) P Edf(^lF) m . 
By Proposition 13.11 letting e — > + on both sides yields 

E^([0,lf) m > (Eexp{-i||Zf||^}) P Eaf([^in m 

In view of (|5.17p . 

liminf — log ( (m\)- Hd ^E\a H ([0, l] p Y 

m— s>oo m V L 



(5.22) 



> 



d(p - 1) log c H + liminf — log ( (m\) Hd(p - 1) E \a H ([S, 1} P ) 



To establish the lower bound for (|5.20p . therefore, it remains to show that 



lim inf lim inf — log - — -ttj-, — -r 

5^0+ m-s>oo m {m]) Hd{ -P^ 1 > 



E 



a H ([8,ir) m >C(H,d,p). 



(5.23) 
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Write 

d H ([0,lf) =a H ([6,l] x [O.ir 1 ) +a H ([0,S] x [O,^" 1 ). 
By triangular inequality, 

l/m 



U[a H ([QAY 



< lE 



a H {[S,l] x [0,1] 



p-l\m 



l/m 



+ ^E 



a H ([0,5] x [0,1] 



l/m 



Given e > 0, 



E 



«f([*,i]x[o,irT 

dx± ■ ■ ■ dx m 



< 



/ ET[p e (W H (s k )-x k )ds 1 ---ds 

J[s,vr k=1 

[ Ef[ Pe (W H (s k )-x k )dsf-ds, 

J\0,l] m f~ 

„ m 

/ ET[p e (W H (s k )-x k )ds 1 ---ds 

Ml]™ t=i 



'[o,i]- fe=1 
da?i • • • dx, n 



p~\ i/p 



} 



dx\ ■ ■ ■ dx r , 



ET[pJW H {s k )-x k )ds 1 ---ds T , 

[0,1]™ ^ 



p~\ (p-l)/p 



E 



i/p 



E 



[af([o,ip 



(p-i)/p 



Letting e — )• + yields 



E 



d ff ([<5,l] x [0,l] p - 1 ) m 



< < E 



Similarly, 



E 



So we have 



a 11 



([0,5] x [0,1] 



p_lym 



< <^ E 



a 



i/p 



E 



a 



H [[o,i] p 



(p-i)/p 



Vp 



E 



a ff ([0,lf) 



(p-1)/p 



E 



a H ([0,l] p ) 



l/mp 



< < E 



^([<UF)' 



a 



l/mp 



+ <^E 



l/mp 



By scaling, 



Thus 



E 



a H ([0,5] p ) m ] = s(p-Hd(p-i))™ E \n H ([0,lf) 



E 



a^([5,lf) m > l-f%-i)// P E(i fl([o il ]P) 



Therefore, (I5.23P follows from fj5. 15f> . 
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To bound the limit in (|5.20p from below, we claim that 



lim — log < (ml) 

m— too ui 



> 



plog |(1 - Hd/p*)-^- Hd l pt \p ¥ )4-(2^y4- J™ (1 + t 2H )~ d/2 e- t dt^. 



(5.24) 



Let n, • • ■ , r p be i.i.d. exponential times independent of B^(t), • • • , B^ (t). Given e > 



E 



of ([o,n] x ••• x [o,r p ])' 



dx\ ■ ■ ■ dx m QP{x\ , • • • , Xjji) 



where 



r°° r /" m 

Qe(xi,- ■ ■ ,x m ) = e _t / dsi • • • ds m ET[ Pe(B H (s k ) - x k 
Jo J\o,t] m t ~> 



'0 U[0,t]- k=l 

Let f(x\, ■ ■ ■ ,x m ) be a rapidly decreasing function on (M. d ) m such that 

/ \f(xi, ■ ■ ■ ,x m )\ p *dxi ■ ■ ■ dx m = 1 . 

J(R d ) m 

By Holder inequality, 

E[af([0,n]x...x[0,r p ])J } 

/ dx\ ■ ■ ■ dx m f (x\, • • ■ , x m )Q e {x\ , • , x m ) 

dx\ • • • dx m f (x\, ■ ■ ■ , x m )H s e (xi , • , x m ) 



(it. 



> 



[o,*] r 



dsi • • • ds m (it, 



where 



fls )e (xi,--- ,X m ) = E JJp^-B^Sfc) -Xfe) S= (Si,-- - ,S m ). 
fc=l 

Consider the Fourier transform 

/(Ai, • • • , A m ) = / dxi ••■ dx m f(xi, ■ ■ ■ ,x m ) exp \i V) A fc • x k \ . 

J(R d ) m I ^ J 

It is easy to see that 

fls, e (Ai, • • • , V> = exp - | £ |A fc | 2 - -Var ( ]T \ k ■ B H \s k )) ■ 

^ k=l ~ k=l ' 

By Parseval identity, 

/ dx\ • • • dx m j *{x\ , • • • , x m )H s e (xi , , Xjji) 

J(K d ) m 

d\\ ■ ■ ■ dX m f(Xx, • • • , A m ) 



(2TT) md J m » 

s m ^ m >. 

x ex p - \ £ i a *i 2 - 2 Var ( £ Afc • 5 " (Sfc) ) 
fe=i fe=i ^ 



31 



Thus, 



{- -I m 1 /p 

E[af([0,r 1 ]x---x[0,r J J)J } 



> 



(2vr 



e t dt I ds\ ■ ■ ■ ds r 

[0,i] m 



d\\ ■ ■ ■ d\ r , 



xf(X lr .. ,A m )expj -|^|A fe | 2 --Vax(^A fe -B H ( Sfe ) 

^ k=l k=l 

We now let e — > + on the both hand sides. Noticing that the left hand side falls into an obvious 
similarity to (I5.7D , 



{- i in -i J 

E[a H ([0 J r 1 ]x---x[0,r p ])J J 

1 



> 



(2vr 



ds± ■ ■ ■ ds m 



[0,*]' 



dAi • • • ciA r , 



f 1 m 1 

x /(Ai,-- - ,A m )exp] -~Vax(J2*k-B H (s k y) 



We now specify the function f{x\, • • • , x m ) as 

m 

f(xi,--- ,x m ) = C m Y[pi(x k ) 



where 



We have 



k=l 



d d(p*-l) 



/" f 1 

/ dAi • • • d\ m f(Xi, ■■■ , Xm) exp <^ - -Var ( V" X k ■ B H (s k 

« f 1 m m 

/ <*7i • • • d 7m exp I - - ^ ~/l - -Var ( ^ 7fcj B^ \s k 



where Bff (t) is an 1-dimensional fractional Brownian motion. 

Let £i, • • • £ m be i.i.d. standard normal random variable independent of Bq (t). Write 

Vk=ik + B^(s k ) k = 1, • 

We have 



, m. 



j lit, 

-Var ( ^ 7 fc% 



k=l 



fc=i fe=i 
And thus by Gaussian integration, 

d7i ■ ■ ■ d 7m exp J -y^Tfe" g Var ( E Tfc-Bofafc)) > 

k=l ~ fe=l ' 

(2vrr/ 2 det{Cov(r ?1 ,-.. ,r? m )}" 1/2 . 
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(5.25) 



with convention that sq = 0. 

Write sq = and assume s± < ■■■ < s m . By Lemma I3~8 



m 

det jCov (771, •• ■ , Vm) } = Var (771) Var 

fc=2 

m 

= {l + Var (Bo(si))} {l + Var (<(s fc )|<( Sl ), • • • ,B^(s fc _i))} 

fc=2 

m 

n{i+(^-^-i) 2H } 



fe=i 



where the last step follow from the computation 

Var (Bg(a k )\Bl*(8i), • • • ,£^-1) 
= Var (b* (s k ) - B^{s k ^)\B^{ Sl ), • • • , B^(s k ^] 
< Var (Btf(s k ) - S^fe-i)) = (s k - s k _i] 

Summarizing our argument since (|5.25p . we obtain 

1/p 



\2H 



{l [«*([(), n] x ••• x [0,r p ])] m } J 



>m!(C(2vr)- d / 2 )' 
= m\{C(2T,r d ' 2 y 



e~ l dt 



1 + t 



dsi ■ ■ ■ ds m TT { 1 + ( 

-| m 

^-^e-'dt 



S k — Sfc-l 



,2H 



-d/2 



I JO 



Equivalently, 



E 



a H ([0, n ] x ••• x [0,r p ])j m > (m\)P(C(27T)- d / 2 ) mp 



'l + t ™yd/2 e -t dt 



pm 



On the other hand, with obvious similarity to f|5. 14[) 



E 



a H ([0,n] x ••• x [0,t p ]) 



< E 



O: 



([0,1]*)] {r(l + m(l-fTd/p*)} . 



(5.26) 



Hence, (15. 24j) follows from (I5.26P and Stirling formula. 

By (I5.16P and (I5.24p . the limit given in (I5.20P is finite. By Lemma [3.91 the large deviation given in 
(EU) holds with 



K(H,d,p) = Hd(p- l)exp 



C(H,d,p)-d(p-l)logC H \ l/H^ (l r , x 

)=c H K(H,d,p). 



□ 
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6 The law of the iterated logarithm 



We will prove Theorem I2.5l in this section. Due to the similarity of arguments, we will only establish 
(|2,18p . By the self-similarity property ([1.121) . the large deviation limit of Theorem 12.31 can be 
rewritten as 

Um(loglogt) -1 logP(a H ([0,t] p ) > \t p ~ Hdl ~P~ 1 } (log log t) Hd ^ p ~ 1 ^ \ 

= -K(H,d,p)\ p */ Hdp (A>0). (6.1) 

Therefore, the upper bound 

limsupt^ p - 1 )- p (loglogt)-^ p " 1 )5^([0,tf) < K(H,d,p)- Hd ^ a.s 

t— >oo 

is a consequence of the standard argument using Borel-Cantelli lemma. 

To show the lower bound, we proceed in several steps. First let N > 1 be a large but fixed number 
and write t n = N n (n = 1, 2, • • • ). Define the d-dimensional process 

ct 



QU(t) = (t + u) H - l ' 2 dB(u) t > 0, 
J 



where B(u) is a standard d- dimensional Brownian motion. Recall that H[0, T] denotes the RKHS 
of {W H (t)} t ^Q^}- Combining Propositions 13.31 and 13.51 we can deduce that {Q^ "(*)}te[o,T] i s n °t m 
H[0,T], T > 0. For that reason, similarly as in Proposition [33J we define the following modifications 
of Q**(t). When H G (0, 1/2), put 



A n t, 0<t<t n 

Qn{t), t>t n , 



where A n = t~ x Q%(t n ). When H G (|, 1), put 

where B 1>n = 3t~ 2 Q^ (t n ) - t~ l Q^ (t n ) and B 2>n = -2t~ 3 Q^ (t n ) + t~ 2 Q^ (t n ). 
Lemma 6.1 For every n > 1, P ({G^ (i)}te[o,t„+i] C H[0, = 1- Furthermore, 

'>*n+l] 



B hn t 2 + B 2 , n t 3 , 0<t<t n 
Qn(t), t>t n , 



S upE||G^||2 , < oo. (6.2) 



Proof: Obviously, it suffices to consider the case d = 1. The first part of the lemma follows by 
the same argument as in Proposition 13.51 For the second part we use Lemma 13.61 with a = t n and 
T = t n+ \. A constant C > below will depend only on H but it will be allowed to be different at 
different places. 
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First consider H £ (0, 1/2), so that m = \H + 1/2] = 1. In this case we get 

7) J 



and 



Therefore, 



and 



/■tn + l ft 

/ / {t- s y H - l/2 G*{ s )ds at 

rtn+i ft n / ft \ 

= C ( (t-s)- H - 1 / 2 (s + u) H - 3 / 2 ds)dB(i 

Jtn JO \Jt n J 



dt. 



E 



tn + l 



(t- s)- H - 1/2 G%{s)ds 

■tn + l ft- 



dt 



= C 
= C 

< c 

< c 



(t-s)- H - x l 2 {s + u) H -V 2 ds) dudt 

*n+l ttn {t _ tn y-2H {u + tn) 2H-l 





" ( I 

1-2H 

n+1 



(t + U) 2 
X-2H ^ 



-du dt 



„dt<C 

t 2 ~ \ t r 



tn+l X~ 2H f tn+1 1 



dt 



tn I t n 



log 



Using bounds (|6.3p - (|6.4p with Lemma [3.61 we get 



(6.3) 



(6.4) 



E ll«n IIh[0A +1 ] < C(*n+1 " tn)tn + C^+l/in)^ \og(t n+1 /t n ) 

< C(N 2 ~ 2H + N l - 2H log N) , 

which proves (|6.2p in the case H £ (0, 1/2). The proof in the case H £ (1/2, 1) follows the same 
line of computations, thus is omitted. □ 

For a simplicity of notation, from now on write M n for H[0,t n +i]. Define the sigma field 

T t = a[{Bx{s),--- ,B p (s)); s < t}. 

To complete the proof of Theorem 12.51 i.e., to establish the lower bound in (|2.18p . it is enough to 
show that for any A < K(H,d,p)^ Hd ^ 1 ^ there is an N, sufficiently large, such that 



"£F{a H ([2t n , Wif) > ACf ^^(loglogWl)^- 1 ^^} 



oo a.s. 



(6.5) 



Indeed, by [9| Corollary 5.29 p. 96], (|6,5p implies that 

hmsup^- 1) - p (loglogt n+1 )-^- 1 )a H ([2t n ,t n+1 ]P) > A a.s. 
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which leads to 

limsup^ d(p - 1) - p (loglogt)-^ p - 1 )d H ([0,t] p ) > A a.s. 

t— >oo 

Letting A — > K(H,d,p)~ Hd ^ p ~ 1 ^ on the right hand side yields the lower bound as claimed. 
Let now e > be fixed and write 



af([2t n ,t n+ i]f)= / ds 1 ---d Sp g e (Wf{s 1 ),---,W p H (s p )) 

J[2t n ,t n+1 ]P 

= [ dsx--- ds p g e {W?(t n + si), ■ ■ ■ , W p H {t n + Sp )) 

= f d Sl - ■■ds p g e {Y l H {s 1 ) + Z? ( Sl ), ■ ■ ■ ,Y p H {s p ) + Z*{s p )), 



[tn,t n +l — tn] P 

where g e (x\, • • • , x p ) is given in (|5.4B and for j = 1, • ■ • ,p, 

ftn+t ftn 

Y 3 H (t) = (tn + t- s^-WdBjia), Zf(t) = {tn + t- sf-^dB^s) . 

Jtn JO 

Consider a symmetric set A C (S>j =1 c|[0,t ri _|_i],R a! | defined by 
A= | (/i,--- ,f p ) € ^ =1 c{[0,t n+1 ],R d }; 

/ tfai---cfa p (/ e (/i(ai),--- ,/ p ( 5p )) > ACf ^(loglogtn+i)^^ 1 )} . 

j[t n ,t n+ i-i n ]p J 

For any (/i, • • • , / p ) 6 <S>j =1 IH n , applying Lemma |3. 71 (h) to the indicator of A leads to 

P{ (Wf + /i, • • • , W* + / p ) e A} > exp { - I Yl 11/11*. } p { ■ ■ ■ , ) € A}, 

if A,-" - ,/ P € H n . 
Notice that 



tn < t < t 



{z H (t); t n <t< t n+1 } = \q* (t); t n <t< i n+ i} = {Gn {t)\ 

{Y H {t); t n <t < t n+1 } = {\V H (t); t n <t<t n+1 } 
and Y (t) and Z H (t) are independent. By Lemma [6.11 

p{ (if + zf , • • • ,y p h + z*) e A\T tn ) 

* ex P { " \ E 11°^ US. } p { (wf . • • • > O e a} 



n+l 



i=i 
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or 



P{af ([2t n ,t n+1 ]P) > A*£f ^(loglogWx)^- 1 )^} 

^ eX P { " \ E lKiHL}P{fif ([*», Wl - *n] P ) > ACf ^^(loglogW!)^ 1 )}. 



i=i 

Letting e — >• + on the both sides yields 

¥{a H ([2t n ,t n+l r) > A<;f ^^(loglogWO^- 1 )!^} 



> exp { - \ £ WGl^}^ {[t ni t n+1 - t n f) > ACf ^^(loglogWi)^ 1 )}. 



By ()6.ip and by an argument similar to the one used for (|5,23p . for A < K(H, d,p) Hd ^ p ^ and any 
small 5 > 0, one can take N sufficiently large so that, for large n, 

F{a H ([t n ,t n+1 -t n f) > XtZ H 1 dip - 1 \loglogt n+1 ) Hd ^} 

>exp{ -(l-5)loglogt n+ i} = (n log N)~ 1+s . 

To establish (|6.5p . therefore, it suffices to show that for any e, 5 > 0, 

1 P 

E E H G »A < eloglogt n+1 } = oo a.s. (6.6) 

n j=l 

Indeed, by Lemma 16,11 can be viewed as a Gaussian sequence taking values in H n . By the 
Gaussian tail estimate, see |31j . p. 59, there is u = u(e) > such that 



1 

v{J2\\<A>^o g t n+1 }<- 

for large n. Then for < 5 < u, we obtain (I6.6p . which yields (|6,5p . The proof is complete. □ 



7 Local times of Gaussian fields 

We begin with mentioning the work of Geman, Horowitz and Rosen (|21j) on the condition for the 
existence and continuity of the local times of the Gaussian fields, see also recent work of Wu and 
Xiao [35]. Let X(t) (t £ be a mean zero Gaussian field taking values in M. d such that there 

is a 7 > such that for any t > and m £ N, 



dsi • • • ds m / dXi ■ ■ ■ d\ m (7.1) 

x ( II l^l 7 ) exp | - -Var ( A fc • X(s fc )) } < oo. 
fe=l fe=i ' 
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Geman, Horowitz and Rosen (Theorem (2.8) in |21j ) proved that the occupation time 

Ht(B) = [ l {X (s)eB} ds BcR d 
J[o,t] 

is absolutely continuous with respect to the Lebesgue measure on M. d . Further, the correspondent 
density function formally written as 



a([0,t],x) = f 5 x (X(s))ds 
J[o,t] 



is jointly continuous in (t, x). For fixed x, the distribution function a([0, t], x) (t G (M + ) p ) generates 
a (random) measure a(A,x) (A C (K + ) p ) on (M + ) p which is called the local time of X(t). 

In this paper, the result of Geman, Horowitz and Rosen is applied to the following four Gaussian 
fields: 

1. The d-dimensional fractional Brownian motion X\(t) = B H (t). 

2. The (i-dimensional Riemann-Liouville process -X^t) = W H (t). 

3. The d(p — l)-dimension Gaussian field 

-X"3(*i, • • • ,t p ) = \ Bx(t\) - B2 (£2), • • • , Bp^tp-i) - Bp(t p )Y 

4. The d(p — l)-dimension Gaussian field 

X 4 (h,--- ,t p ) = (w 1 H (t 1 )-Wi I &),■■■ ,W^i(t P -i)-T<(t P 

Theorem 7.1 Under Hd < I, Xi(t) and X%(t) satisfy the condition l\7.1\) : under Hd < p* , -X^t) 
and X&(t) satisfy the condition \7.1}) . Consequently, X\, X2, X3 and X4 have continuous (jointly 
in time and space variables) local times. 

Proof: Due to similarity we only verify (|7.ip for X3, which becomes 

/ dsi---ds TO / dXi ■ ■ ■ d\ m exp \ - -Var ( V" A fc • X(s k ) ) \ TT (Afcl 7 < 00, (7.2) 

where we use the notation 

Sfc = (sfc,i) • • • > Sk,p) an d A& = (\k,u • • • , \k,p-i) ■ 

Notice that 

m p m 

Var (JT Ajt • X(s fc )) = ]T Var ( £(A fcJ - A w _i) • B H (s k ^ 

k=l j=l k=l 

with the convention A/^o = Afc )P = 0. By suitable substitution and using the bound 

v 

|Afc| < C J|max{l, |A fcj - - Xk,j-i\}, 

3=1 
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we have 



/ d\i ••• cZA m exp <^ - -Var ( V" ~\ k ■ X(s k )) \ TT |A fc |' 

j (Rd(P -i ))m [2 yj 11 

r p 

<C / dX 1 ---dX p -x[\H j (X j ) 

J( R md)p-l 11 



where 



flj-C^) = ( II max ^' I'M 7 }) ex P \ - 2 Var ( E M • B ^Ofc,i)) f 

k=l ' fe=l ' 

for A, = (Aij, • • • , Xmj) (1 < i < p - 1) and A p = -(AH h A p _i). 

Write ^ 

n^(A,)=n n ^(a*) 1 ^. 

3=1 j=l l<k^j<p 

By Holder inequality 

f dX 1 ---dX p „ 1 T\H j (X j )< TT { / dAi-.-dVi TT ^fe( A fc) p 

' / ( Kmd ) p - 1 7=1 t=l I (K md ) p_1 Kfcii< B 



When j = p, 

/ dAi---dA p _i IT ^fe(Afc) p * = IT / ^(A) P *^A. 

• / ( Rmd ) p - 1 i<fc< P k=i Jvnd 

As for 1 < j < p — 1, recall that X p = — (Ai + • • • + A p _i). By translation invariance, 

/ H p (X p f dXj = f H p {Xf'dX. 

By Fubini theorem, for fixed j, 

[ dX^-dXp-x TT H k (Xkf = TT I H k (XfdX. 

■ /(RW)p " 1 i<k£<p i<k£<P jRmd 

Summarize our argument, the left hand of (|7.2p is bounded by 



k=i 
l/p* \ P 



{,. r „ m 

/ dsi---ds m / dAi---dA m (rfmax{l,|A fc | p ^}) 
J[0,t] m lJ(R d ) m K Z 

xexp -P-V ar ^X k -B H (s k )) 

^ k=i ' - 

Hence all we need is to find 7 > such that 



ds\ • ■ ■ ds r 

[o,t] m 



/ dAi---dA OT (jI|Afc| 7 )expi -|-Var(^A fc -i? // ( Sfc ) 
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for all m G N. Further separating variable and substituting variable, the above is reduced to 

d/ P * 



I 

J[o,t]* 



dsi ■ ■ ■ ds r 



n lib f 1 

/ dAi---dA m (n|A fc r)expi _-Var(^A fc ^( Sfc 
-' Rm k=i ' ^ fe=i 



< oo. (7.4) 



By (|48|) . for any = s < si < ■ ■ ■ < s k , 

Var [Bg(8 k ) - B^(a fc _i)|B^(ai), • • • , B^Wi)) 

> 7^(s k - s k ^) 2H = ^Var (B*(8 k ) - B?(a fc _i)). 

This property is generalized into the notion known as local non-determinism. By Lemma 2.3 in 
Berman [8], there is constant c m > such that for any Ai, • • • , A m G K and any si < • • • < s m 

Var A fc (5 % fc ) - 5^(a fc _!))) > c m ^(s fc - s k ^) 2H X 2 k . 



k=l 



k=l 



Consequently, with notation Aq = 0, 



/ dAi ••■ dA m ( J] |A fc r) exp - -Var A fc Bf (a fc )) 
•^ Mm fe=i ' fc=i ^ 

« m 

= / dAi • • • dA m ( TT |A fc - A fc _i| 7 ) 
Jw> \ k=i J 

( 1 m 1 

xexp --Vai(£A fc (B^( aj 0-B?(* fc -i))) 

^ fe=i ^ 

,, m ^ m 

< / dAi • • • dA m ( TT |A fc - Afc_i| 7 ) exp \ - c m V^Sfc - s k ^i) 2H X 2 k \ 
J ® m fe=i 7 L t^i J 

Using triangle inequality (for which we take j < 1) 

mm m 

ni A fc-^-ir < n (iA fc i 7 +iA fc _ii 



fe=i fc=i 

where 8a = 0, 7 or 27. Notice that 



n < n a v \ x k\) 6 ^ < Ha v iA fc i) 27 . 



fe=i fc=i *!=i 

Notice the number of the terms in the previous summation is at most 2 m . Thus, 

m m 

U\ x k- x k-ir<2 m Yi(iv\x k \) 2 ^. 

k=l k=l 

In this way, the problem is reduced to finding 7 > such that 



'[o,t]: 



dsi ■ ■ ■ ds r 



n lib f lib 

J d\i ■ ■ ■ d\ m ( If |Afc| 7 ^) exp < -c m ^( 
-* Rm k=i ^ k=i 



{s k - s k ^) 2H X 2 k 



d/ P * 



<oo. (7.5) 
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Observe that 

-c m (s t -sn) 2fl A 2 ^ 



d\i ■ ■ ■ d\ m ( 11 \\ k \A exp l-cm^k- s k -!) 2H \i = II / l A l 7e " 

fe=l k=l > k=l J -°° 

-oo J 



k=l 

Therefore, we need to choose 7 > such that 



„ m 

/ cfei ■ ■ ■ cfe m TT ( Sfc - afc.xJ-C 1 "^)^ ' < 00. 

This is always possible because if<i < p*, so that (1 + 7)i2"d < p* for some 7 > 0. The proof is 
complete. □ 



8 Appendix 

Lemma Al Let {B H (t)}teR be a standard fractional Brownian motion given by 

B H {t) = c H f ((t - sf^ 2 - (s)^ 2 ) dB(s), (Al) 
J— 00 

where {B(t)}t£R is a standard Brownian motion. Then 

c H = V2H2 H B(l - H,H + l/2y 1/2 , (A2) 
where B(a, b) = x a ~ l (l — x)^ 1 dx is the usual beta function. 

Proof. Since Y&r(B H (1)) = 1 we get 

-1/2 



CH 

Put 



{[ { {i+x)H ~ i/2 - xH ~ i/2 ) 2 dx+ ^Y ■ (A3) 

1 = J [(l + xf'^-x 11 - 1 / 2 ) dx. 



Then 



/•oo , . 

/= lim / ((l + xf-W-x 11 - 1 ' 2 ) e-^dx 
= lim ((e" + l)^ 2H T{2H)-e^^ 2H 1 {2H^)-2 f°° (1 + xf^x^^e'^ dx\ 

M^0+ [ Jo J 

= — — + lim !.2e^ 2 fi- 2H T(2H) - 2 f°° (1 + x^'^x^^e'^ dx\ 
2H ^0+ I Jo J 

= " ^7 + & {*"V"iW - A e ,/ 2r(ff + 1) } , 
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where 7(2, x ) and K v (z) are the incomplete gamma function and modified Bessel function of the 
second kind, respectively. The third equality uses the facts that e ll iJ,~ 2H j(2H,fi) = + o(l), and 
that (e** + l)^ 2H = 2e^l 2 ^ 2H + o(l) for fl" < 1, as /i — > 0. The forth equality applies formula 
3.3838 in [22]. 

Using the duplication formula T(2H) = ¥^J-T{H)T(H + i) (see [22, formula 8.3351]), we get 
Since 

fl -2H 2 2H T{H)= f°° x H-i e -^dx and K u (z) = -(-Y r t^e-'-idt 

Jo 2 2 y 

(see |22| formula 3.4326]), we obtain 

1 1 „/„ 1 



00 

J= — — + -^=T(H + -) lim / x H - 1 e-'r x (l-e~&)dx 
2H v 2%^o+io 



00 

x H_1 (l -e~^)d:c 



2# 

- 1 1 r(i-ff)r(F + i) 

-"2tf + ^F4^ • (A4) 

r(x)r(y) 

Combining (IA4I) with (IA3I) and using the well-known formula B(x,y) = — (see, e.g., 1221 

T{x + y) 

formula 8.3841]), we get (|A"2~j) . □ 
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